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Abstract 

It is shown by means of a simple analysis that the linearized system of 
transport equations for a relativistic, single component ideal gas at rest 
obeys the antecedence principle, which is often referred to as causality 
principle. This task is accomplished by examining the roots of the dis- 
persion relation for such a system. This result is important for recent 
experiments performed in relativistic heavy ion colliders, since it suggests 
that the Israel- Stewart like formalisms may be unnecessary in order to 
describe relativistic fluids. 



1 Introduction 

In 1940 C. Eckart published three papers entitled The thermodynamics of ir- 
reversible processes [T], the third one addressing the problem of a relativistic 
simple (single component ideal gas) fluid. In that paper, Eckart proceeded fol- 
lowing the basic ideas of classical irreversible thermodynamics [2j, except for the 
fact that he introduced relativistic terms in the energy-momentum tensor. As 
part of his phenomenological approach, he proposed constitutive equations with 
relativistic corrections. Since Eckart's theory apparently leads to results that vi- 
olate causality and involves undesirable unstable modes [3], it has been patched 
up in several ways using formalisms introduced by Israel and coworkers [1] [S] 
^ and sometimes using extended irreversible thermodynamics ^ ^ . Recently, 
it has been shown that the unphysical behavior of the unstable modes is due 
to the coupling between heat an acceleration proposed by Eckart [5] . Indeed, it 
has been shown that such a relation is not sustained by kinetic theory [lOj. 

The so-called causality problem of heat conduction, which should be more 
precisely stated as antecedence problem, remains still a controversial issue which 
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suggests the need of extended theories. We wish to point out that eventhough 
the term causahty is the most favoured, the problem of faster that hght propa- 
gation of fluctuations is not striclty a cause-effect issue but rather the prediction 
of an unphysical behavior concerning arrival times. However, as is shown in the 
following sections, relativistic classical linear irreversible thermodynamics, as 
obtained from relativistic kinetic theory features both stability for the equilib- 
rium state [10| and satisfaction of the antecedence principle. 

To accomplish this task we divide the rest of the paper as follows. In Sect. 
1.2 we recall the Navier-Stokes equations for a simple relativistic fluid [TT] and 
introduce the appropriate constitutive equation for the heat flux. The linearized 
set of transport equations is thoroughly analyzed in Sect. 1.3. Conclusions and 
final remarks are included in Sect. 1.4. 

2 Transport equations for the relativistic single 
component ideal gas 

The starting point are the balance equations for a relativistic fluid which are 
obtained from the conservation of the particle density flow 

= nv!" (1) 

and the energy- momentum tensor which, following Eckart[T] reads 

T/^ = ^w^u^ + phi^ + TT^ + \q''u^ + ^u^'q^ (2) 

In Eqs. ([I} and ©, n is the particle number density, u'^ the hydrodynamic 
velocity four vector, c the speed of light, p the hydrostatic pressure and = 
(5^ 4-u''Mi//c^ the spatial projector. The internal energy per particle, e , includes 
the rest energy since it is given by |14j 

£ = mc^ ( 3z + I ^ mc^ -f \kT + ... (3) 

where z = is the relativistic parameter and /C„ (i) are the modified Bessel 
function of the second kind. The dissipative fluxes are the Navier tensor tt^ and 
the heat flux . The conservation equations iV.^ = and T^^ ~ for the 
quantities defined above yield the Navier-Stokes equations for the relativistic 
simple fiuid namely, 

h + n6 = (4) 

C^^ P \ ■ f ^'^ P r\ in u 

+ {qt^^u^ + q^'u^^^ + Oq^ + u'^g^,^) = (5) 
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nC„T + (^^^ e + + + l^t-^g, = o (6) 

where kt is the isothermal compressibihty, j3 the thermal expansion coefficient 
and Cn the heat capacity at constant particle density. As it will be remarked 
below, the tensor tt^ can be further decomposed in a traceless symmetric part 
and the trace multiplied by the spatial projector. 

In order to close the system of equations, constitutive relations for the heat 
flux and Navier tensor must be introduced. The equation for the heat flux has 
been recently established by means relativistic kinetic theory and reads 

c{ = ~Lt'— + u'^ (7) 

where Lt and L„ are transport coefEcients|10|. A detailed discussion on this equations 
can be found elsewhere [10|. The equations for ■k'^ in Eq. ((2} are well-known namely, 

= -2r)a^, (8) 

tr (tt) = -Sy ■ u (9) 

where Tr^t' is the symmetric and traceless part of Tr^f, tr (tv) its trace and a^i, is the 
symmetric and taceless part of the velocity gradient. The transport coefficients in Eqs. 
(O and ([9]) are the shear and bulk viscosities respectively. 



3 Linearized relativistic hydrodynamics 

In order to linearize the set of equations we consider n = jiq + 6n, T = 

Tq + ST and u'^ = du'^ where naught subscripts denote equilibrium quantities 
and the S prefix indicates small perturbations around it. With this hypothesis, 
the linearized transport equations for a simple, relativistic fluid in the absence 
of external fields are 

Sh + noSe = (10) 



— {noEo + Po) 5uu H 5n H 5T ^ 

-(Se^,-27j{S<j^;)^^-^ST,-^^Sn,, =0 (11) 

nCnSt + (^^^ 60 - {LtST''' + LnSrv''),!^ = (12) 

where we have defined 9 = u^^. It is important to point out that the transport 
coefficients in general depend on the state variables. However, since they only 
appear as factors of derivatives of the corresponding fluctuations, considering 
fluctuations on them would induce higher order terms, which are neglected in 
the linear approximation. 

It is crucial at this point to make the following observation. The so-called 
causality violation of the transport equations to first order in the gradients, given 
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by linear irreversible thermodynamics, can be easily spotted by observing that, 
considering Uq = and linearizing, Eq. ([6]) leads to a parabolic equation for T. 
This clearly admits arbitrary propagation speeds for the corresponding signals. 
However, the hypothesis of a fluid at rest or the fact that calculations can be 
performed in the comoving frame should not be translated into a vanishing 
hydrodynamic velocity, but in Uq = as considered above. That is, Su'^ should 
not vanish even for the fluid at rest or in the comoving frame; only the mean or 
equilibrium velocity can be zero. This fact has already been pointed out in the 
analysis of the linearized relativistic Euler regime [H]. 

The analysis of the dynamics given the system of equations (|10m2p can be 
found in detail in Section 4 of Ref. [TT] where we discussed the modifications 
to the Rayleigh-Brillouin spectrum. Here we only quote the results needed in 
order to address the problem at hand namely, the causality of the system. We 
start by calculating the divergence of Eq. (fTTj). The transverse mode is then 
uncoupled from the system and a set of three scalar differential equations for Sn, 
SO and ST is obtained. A Fourier-Laplace transform is then performed, leading 
to a system of algebraic equations depending on the time and space variables, 
s and q respectively, whose associated determinant reads 
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(13) 



nocn ^ noCrtKT nocn ^ 

where for convenience we have introduced the following notation[ll) 

Po = (n-o£o + Po) (14) 

A^C + 477/3 (15) 
The dispersion relation is thus given by 



s 



3 + d2S^q^ + s [d^q^ + diq^) + d^q^ ^ (16) 



where the coefficients d2 to have been specified in an earlier workjil]. The 
physical interpretation of the three roots of Eq. ([16]) is well known. The dy- 
namics of the perturbations in the fluid are characterized by a strictly dissipa- 
tive component which decays in time depending on the value of the real root 
while the other wave-like component propagates at a speed given by the imagi- 
nary parts of the conjugate roots, damped by a coefficient which depends on a 
Stokes-Kirchhoff like factor. Moreover, a plot of the dynamic structure factor 
as a function of s for a fixed g, will feature three peaks. In this work, we are 
interested in the location of the symmetric Brillouin peaks[13j, which are given 
by the imaginary part of the conjugate roots, that is w = ±\/dAq. Thus, in this 
case 

^ = i^/^? (17) 
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such that, the distance between the peaks, i. e. the speed of propagation of the 
wave-Uke component of the fluctuations, and the origin is bounded by 
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KT [noeo + po) 



(18) 



Notice that, in the non-relativistic case, the fluctuations propagate at the speed 
of sound, i. e. 

„2 _ 7 



KTPO 

As an example, for an ideal gas 7 = 5/3 and kt — such that 

, 5 kT 
6 m 



(19) 



(20) 



which is clearly unbounded and can be increasingly large for high temperatures. 
However, the speed of propagation in the relativistic calculation reads 
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KT (noEo +Po) 



(21) 



Using the expression for the internal energy density given by Eq. ^ we now 
obtain 



3z- 



(22) 
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As can be seen in Fig. 1, the expression in curly brackets never exceeds the 
unity which finally shows that the propagation speed for signals in the transport 
equations for the relativistic fluid does not violate the antecedence principle. 



4 Final remarks 

The transport equations derived from Eckart's modified theory show, in its linear 
version, no problems regarding stability and causality features. This fact implies 
that there is no real motivation to introduce the so-called second order theories, 
which introduce non-fundamental adjustable parameters. A simple first order 
formalism is desirable to describe the fluids which are formed in RHIC type 
experiments [ 1 5J . 

In the non-linear case the relativistic Navier-Stokes equations here employed 
are, by far, more complex. In this context, very little can be said regarding the 
problems of stability and causality with the techniques included in this paper. 
It is desirable to perform further work in this direction. 
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Figure 1: The ratio of the speed of propagation to the speed of hght, as a 
function of the relativistic parameter z. 
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